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FOREWORD 


The value of microwave scafferomelers and radiomefers as remote sea wind 
sensors has been independently demonstrated by a number of investigators. However, 
near-simultaneous observations by a composite- radiometer and scatterometer 
(RADSCAT) instrument have been judged to have value in making better estimates 
of the surface winds beyond the improvement provided by two independent measurements. 
To demonstrate this potential a joint effort between New York University, General 
Electric Space Division , the University of Kansas and NASA Langley Research was 
undertaken through the Advonced Applications Flight Experiment program of NASA. 

This document reports the investigations performed by the University of Kansas during 
the first year of this joint program. 

Specifically, this report was prepared by the Remote Sensing Laboratory of the 
University of Kansas Center for Research, Inc., under Contract NAS 1-10048. 
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I . INTRODUCTION 


The microwave emission characteristic of the sea has been measured by several 
investigators [1,2, 3, 4, 5]. These investigators have compared their observations 
with predictions from geometric optics theory [6] which uses a single surface mode! 
and found some but not satisfactory agreement between predictions and measure- 
ments. The wind dependence of the geometric optics approach was based on measured 
rms sea slope data presented by Cox and Munk [7]. However, the theory failed 
to predict the observed emission characteristics near nadir and fitted only loosely 
for nadir angles between 30 and 70 degrees. The failure of the geometric optics 
model to account for wind dependence at nadir was first reported by Nordberg, 
et al. [1], and verified by Hollinger [5]. 

In view of the above deficiencies, an investigation is necessary to seek a 
more adequate model for microwave emissions from the sea. The emphasis in this 
investigation is oriented towards using a composite surface model which better 
reflects the roughness characteristic of the sea. Since several lengthy numerical 
integrations are required to yield the emissivity, the more adequate model must 
not be so complicated as to make numerical calculations prohibitive. With this 
perspective, a non-coherent scattering theory of the type described by Semyonov [8] 
is extended to yield the bistatic scattering coefficient. Since an acceptable scat- 
tering coefficient for predicting the microwave emission characteristics must also 
be acceptable for predicting backscattering, the latter case is examined to provide 
a cross check on the model. 

To compare with experimental observations, an isotropic surface character- 
istic, although not realistic for the ocean surface, is assumed. A justification for 
this assumption is based on the observed directional insensitivity of emissions from 
the sea [9]. The two-scale rough surface model is also assumed to have Gaussian 
surface height distribution and Gaussian surface correlation for both scales. The 
dielectric constant needed in the calculations is based on the data reported by 
Saxton and Lane 1 10]. 

The wind dependence of the surface parameters in the composite model Is 
introduced in accordance with rms slope data measured by Cox and Munk for the 
large undulations and Sutherland's [11] results for the small irregularities. Details 


for the theory and the choice of surface parameters are given in sections li and III 
respectively. Comparisons of the computed brightness temperature and backscattering 
characteristics at two different wind speeds are made with both measured data and 
the predictions of a single surface model. The results are presented in section IV. 
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II . A SURFACE BRIGHTNESS TEMPERATURE THEORY 


11.1 Formulation of the Problem 


The basic theory of surface brightness temperature was developed by Peake 
[12]. The relationship among the surface emissivity, the surface temperature and 
the brightness temperature is as follows: 

T BjC0J=ej (0 ) Tg j = k or -ir 

whereTejO) is the brightness temperature;€j( 0 jthe emissivity; Tg the surface temper- 
ature; k denotes horizontal polarization and ir vertical polarization. Note that 
the azimuthal angle <f> needed together with the nadir angle 6 to specify a direction 
has been chosen to be zero without loss of generality. 

The connection between the emissivity and the differential scattering coef- 
ficient of the surface, 

€j(6) = I yj (0,0s,1>s) d0sd‘l>s (1) 

where 0 s,<^> 3 are angles defining the direction of scattering corresponding to a wave 
incident at an angle 6 . 

The basic formulation of the problem indicated above shows that the differen- 
tial scattering coefficient is the quantity that defines the angular characteristics of 
the brightness temperature of a given surface. Consequently, different brightness 
temperature theories are also distinguished by the different models assumed for 
the differential scattering coefficient. 

Under the non-coherent assumption 65 , <l> 5 )can be shown [ 8 ] to consist 
of two terms, I .e . 

yjC0,0s,<^>s) = + < yj\o,ds,<Ps) > (2) 

where denotes the main contribution by the large undulations; 

denotes the differential scattering coefficient of the small Irregularities averaged over 
the distribution of the surface normals of the large undulations. Detailed derivations 
for and^3j(6^9j,<[’i))are given in sections 11.2 and 11.3 respectively. 
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Since the backscattering cross-section is a special case of the differential 
scattering coefficients, it can be obtained from the knov/n differential scattering 
coefficients, (see Figure 1 ) i.e. 


o-.-(e) 

or equivalently 


COS 0 




0s = 0 

«i‘s=Tr 


C&) = Ce) 1- < c6 ) > 


(3) 


<^3° (8) = cos 6 • a'/ (6, 0s >5) 6 ^=e 

^ *^s = Tr 

<^s^iC6)}=cosQ • 

4*s=«Tr 

Detailed derivations of given in section II. 5. 

II. 2 Derivation of 


(3a) 

(3b) 


To derive ^(6,0 we may begin with the vector scattered field due to a 
plane wave incident upon an undulating surface to which the tangent plane approx- 
imation is assumed applicable. Such a field expression is, in general, rather 
complicated. To simplify results, the stationary phase technique will be employed. 

An expression for not indicating explicitly the effect of the small irregularities 

will be derived first. This expression is then clarified to reflect the small structure 
effects by computing the explicit forms of the modified Fresnel reflection coefficients. 
As pointed out by Semyonov, such a computation may be performed for the more 
general finitely conducting surface in accordance with Rice's paper [13]. 


II. 2. 1 The scattered field 

The far zone scattered field in the direction 2 J, 2 due to a plane wave polarized 
along ^ Impinging upon an undulating surface Z (’f.y) can be shown to be [1 4] 


A S 




E e 




ds 


(4) 
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where a time factor of the form e ^ has been suppressed; K“ — atTr — > R is the 
distance from the origin to the field point, k is the wave number; <.Rh),<Rv) are 
the modified Fresnel reflection for horizontal and vertical polarization respectively; 

is the magnitude of incident electric field; 'S the unit propagation 
vector of the incident field, and 7l is the normal to the surface 


The set of orthogonal unit vectors serving as the local coordinates 

for evaluating the local field on surface is illustrated in Figure 2. The unit vectors, 
"t and d relate to n and /n,as follows 

/VVS 

■fc = C^i X TT ) /) -n, X n 1 

AAA / • AA<A | 

d = 'a , X t 

where may be written in terms of the partial derivatives, Z-x, of Z ( “X, ^ ) 

along x and y axes as follows 

C I 6 , k ) are the unit coordinate vectors. 

*VH AAA 

To simplify (4) by the stationary phase approximation, let 
j ^ be the vector components of ^ . Then the phase factor in (4) is 

I'l = J,’' + ?-j ^ (6) 

The stationary phase assumption requires 

^ ♦ 

It foil ows that 

H, = - 

2 3 = - fj/i, P) 

, ^ l/iii. . 

The signTficanf result indicated by (7) is that all surface slopes in (4) may 
be written in terms of the incident and scattered propagation vectors. Consequently, 
the integrand in (4) except for the phase factor, expf- jk ^ • r T can be moved 
outside of the integral sign, i.e. 
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with 



= EoE: nz ^ ( 


- B i +nj 

. XCC 

D) 

X r jk(^a'2l'3 

• r 

aw 

ds 








(8) 

A = 

0+ <Rk» c 

<x • 

AAA. 

-t ) ( ^ 

>vn '■ AW 

X ± ) 

An ' 




B = 

A^A 

(i- <Rv»( 

a • 

•d ) C71 

Av^ >v\ 

■71.) 

AW 

-t 

xvv> 


(9) 

C = 

ywv\ 

c '+ < Rv> ) ( 

CC • 

AAA 

cJ ) ( n 

J AAA 

X -t ') 

AAA ^ 




D = 

Ci-<R,^>) ( 

(K ' 

-t ) C'n • 

AAA a aaa 

n, ) 

t 

AV« 




The local coordinate vectors, two other vector products in (9) and the differential 
surface element in (8) may all be written in terms of the propagation vectors, 
t = C -n, n .2 ) / 1 n. X ria I 

d ^ .( • 21* ) Tl- - Tka 


nx ± =. 


• -n , = 


1 ^ ri.* I 
- I 21a - n ■ I 


( 10 ) 


ds = ( liia-ri 

To express the polarizd 


l/^z ) d% 


3tion states of both the incident and the scattered j 
field. It Is'jconvenient to introduce a set of orthogonalj unit vectors for | 

the Incident field and another set for the scattered field (see Figure 1 ). In 

view of Figure 1 , explicit expressions for these unit vectors may be written 
- 'll, = - sin 6 i + cos 6 
“V, = COS0 •+ 51.71 9 J< 

k, = j 

= sin 9s 4>s i, + sin 6 ^ ^ 

= - COS 0s cosct.s;i,- cos ©s sin4>s 
K, = sin >i - cos 4’c J 

For horizontally polarized incident wave, fhe polarized and cross 

polarized scattered fields are given respectively by 
^kK" ^*-2 ' 

= EoK [kz AK-BK)t7l2xCC,,^Dh>}] f e ^ ds 

E Kv~ ifa ■ EK(! 


( 11 ) 


( 12 ) 


= E.£to • 2,^ " S'')) I 


(13) 


ds 
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where 


Ak= -i-fi c 

Bk = i-f I ( I- <Rv>U 


Ck = -ifi(it<R„» (h.. 

fik= \i\ ( <Rk» Ci;. -no 

lfl = 


1 23.2 - 


2| 21' * 

Similarly, for vertically polarized incident wave polarized 

and depolarized scattered fields are 

^ vk = s r jh(22,-n..>-X 

= E.KCjb ■ 

ElrV ~ 2^2 ■ ^ 


cis 


V, S 


where 

A,= /f(({ + <R^» 

Bk= IfKi- <Rv» •!}.= ) 

C».= -|fH 1+ < Rv>) (^iO ■ 3i^)^ "■■' 

iv=-|flCi-<RK»Ck.--lO(2>'^n^) 

The above field expressions may be further simplified using the vector 
identities 

• oa 2*( M + 'n 2 X Nl ) = ir, • M - ka • N 




(14) 


dS 


Thus, 


'1^2x(M + n2xN3 =-KE'M-'^ 2 *N 

•^V^A ' ✓'vvs /VV^ /VVv AAA 

Ehh = C.^-<Rk> be + <Rv>df ) I 


Ekv= CoC ■«•)■ ci c ) I 

Evk = Co(-<RK>df +<R^> be) I 
E\r\r'CoC<.RK>dcL+ 
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where 


-ik R / ^ 2 V 

I = da 

%= sen 0s cos4>s - sin 0 
tfj = sin 05 sin 
fjr - cos 6 + cos 6s 

b = ( If' * 23-x ^ = sin. 0 cos 0^ + cos 0 sin 0s cos <|>s 
^ sin. 6 cos 9 i cos - cos 0 sin 65 

^ C i' *I!;0 = sin 0s sin<<>^ 

■f - Sin 0 sin<#>s 

= / 2;2 - 2i'/^ /a = I + cos 0 cos 0s - sin 0 sin 0j cos 

Cc'tf') = CbSd") 

Note that by expressing n^In terms of the orthogonal set (>ii ti.' ) 
it can be shown that 1.2;* ’’ 2:.' I* = b* -r d* . Similarly, by expressing n, in 
terms of the orthogonal set (nj , Vi , JK,.), it also follows that 1 22;i '' Zl* 1** c*+ -J- * . 


11.2.2 The differential scattering coefficients 

The differential scattering coefficients related to the scattered fields computed 
in the previous section [12,15] are of the form 

where the symbol * denotes complex conjugate, and ^ the ensemble average. 
The subscripts j , -i denote the polarization states of the incident and the scattered 
fields respectively. Ao is the illuminated area. 

Upon substituting (15) into (16), it follows that 

i- C0.0s.</>s) 

= r KRh>l%^l<Rv>l'd^ -1 

■^AoCose?-/ I Le^js . J / 


■n- A,cos0 k + 

o 

^v(0, 0s,</>s) = (0,05.<^s) + i>'uv(0. 05 ,<^s) 

= k" < ( KRv>l%^-f HRK>l^d^ 




'Tr Aocos 0 


b=-vd‘ 


-J <m'> 
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For on isotropically rough surface with Gaussian height distribution, III ^ 


simplifies to 

< IH^> = 2 tt AoJ j,(k e 


( 18 ) 


where J„( ) is the zero order Bessel function of the first kind; and J>(5)are the 
variance and the autocorrelation coefficient of the surface respectively. 

An approximate solution for (18) consistent with the stationary phase approx- 
imation is 




z 'rr Ao expf- f'x fn ] 


(19) 


where “W. = [o'^j f (o) | ] ^ is the rms slope of the surface . 


11.3 Derivation of 




The scattered field due only to the small irregularities has been derived by 
many investigators [1 3,16,17]. As indicated in the previous section, with the ' 
scattered field expression known the scattering coefficient can be computed. To 
account for the interaction between the small irregularities and the large undulations, 
the expression for the scattering coefficient of the small irregularities is then averaged 
with respect to the slope distribution of the large undulations. The resulting expression 
is the desired scattering coefficient, 

11.3.1 Differential scattering coefficients 

The far zone scattered field of i-polarization along the direction defined by 
the angles 0s and 4s due to a plane wave of j-polarization with unit amplitude 
impinging upon an irregular surface s(<, 3 )aIong the direction defined by the angles 0 and 
4'’ has been derived by using the method of small perturbation [17]. The ensemble 
average of the magnitude square of the scattered field, <| Eji(6,4.' 6s/4s ) ^ 

can be shown to be [17] 


< (Ey^( 9, 4,' ds',^s^r) - 


cos^O cos^ &s 


( 20 ) 
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where 


j.i,. = incident and scattered polarizations respectively, either horizontal 

or vertical polarizations, 

A = illuminated area, 

R = distance from the field point to the surface, 

(€c- n cos ^ ) 


Mkv “ 

MvK = 
Mv V 


(cos 0'+ -stride; ) 

( e >■ - I ) sen Jer- sLn^0i 

Ceos e'i-J^r-sLn^Q')(e^rCos 6^ + - sLn^ 6 i ) 

- ) sin J^r-sLn^ &' 


( 21 ) 


( € rCOS0'+Jer -SLn^a' )( COS 65 + /6>--Sin®0s ) 

_ ~ 1) (gos i<Ps~4>)hr-slriQ’ si-n^e^ - 6r si7i0g sin0 ' ] 

(CrCOS 6'+y6^-si-n^0 ' )(^-rCOS 0s 7cr- sLti^ ) 


6y = complex relative dielectric constant, 
'W’Cp.'Jl" surface roughness spectral density, 
p = k (sen 6s cosC<^s stn0 ') 


^ = k sen 9s si-n ) 


Note thpt the surface roughness spectral density w(p,q), is related to its correlation 
function R (x,y) by 


R C^, y ) = 




dP* + jf'i- 

w ip,%) e 


d p d<J. 


For an Isotropically rough surface. It reduces to 

R CD = ^ pw (i) J,Ctr) tdt 

and ^ 

-wrc-fc) = ~rr ROrj J^Ctr) rdr 

-'O 

For a Gaussian spectral density, it Is expressed in the form 

w (t) (22) 

where 

■t = 

^ 1 , ^ = standard deviation and correlation length of the small irregularities 
respectively . 
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From (16) and the relation 

y/ ( e ' cjy' 0>,'6s,0s 3 + Jf'jvC 0>,’ 0s .^5 ) 


we get 

= 4 k^q"i"cos0We;(|Mjh|%|M^-vrjexp[-(#) } 

■j = k or V 

11.3.2 Averaging Procedure 

To account for the tilting effect of the small irregularities by the large 
undulations it Is necessary to average C^, 4>'. Qs.'Ps) with respect to the slope 
distribution of the large undulations [8]. That is 

<jr/c0.e„4>5)> = f"’ f” ( a>; e ^,^' > p c 2, , 2 ^ 3 / 1 + z| + z/ dz^dz^ 

^ 00 00 


where 0 Is the incident angle and &s cind are the scattering angles. 

To evaluate the above integral, connecting relations between the surface 
slopes Zx.Zj and the local angles, 0',^^ 0s>‘^’s °>'e needed. To find these 

relations let us first express Zx ond in terms of the azimuth and the elevation 

angles, 4k. and 0,^ r which represent the tilting effect 

= COS tan Bn 
= sin^,^ "tan 6n 
From (25) It follows that 

2^+2', = 5^^ 8" 

0 ' « t On , T-n ) ' 

The local angles (0,4*, ®s,‘4’5 ) may now be related to the azimuth and the elevati 
angles 4k. and 9^ by connecting relations derived below. 


In Figure 3, the two sets of coordinates (x,y,z) and (x', y', 
in terms of the angles 0 n and 4 *tv as follows: 


; 

-X 


cos 0n COS 4^ 

COS 6„ sin 4,1 

sin 6n 

✓ 

3 ' 

- 

- sen ^•n 

cos 4-n 

0 


i’ 


-St;n0„ cos 4n 

- sin sin 4n 

cos 0n 






«< 



z') are related 


7f 

% 


(27) 


Hence, for a scattered field point, P, located at a distance, R, from the origin 
the coordinates of P may be expressed either in terms of the angles G 5 and ^5 
the angles 65 and <^5 , 

-x’ = R sin 6s cos <^5 

3 - = R sen 6s 

2' = R cos 6s 

= R sLn 6s cos 4 s 

^ = R si.n 0 s si.u 4 s 

2 = R COS 6s 


f 

or 


(28a) 


(28b) 


Substituting equation (28) into (27) we obtain the connecting equations for the sets 
of angles as follows: 


sin 0s cos 4 s 


y N. 

cos 071 cos 4 rt cos0n sin 4 -n sin 0n 


f > 

sinflsCos 4 s 

siTi 6s sin 4 s 

- 

- sin. 4 ^ cos 4.^ ° 


sin 0s sin 4 s 

cos 0s 


- sin 0.T, cos 4-71 -sin07iSin4n 005671 

s. ^ 


cos 6 s 


If we take the angles with a prime to be local scattering angles, we see from 
(29) that the local scattering angles may be expressed in terms of the scattering angles 
4 s tilting angles 6 ^ , , i.e. 

COS 0s = 0^ COS0. - SI.U 0-n scnGj cos (<p ~ cl>^) 

( 30 ) 

sc-n = sen 0 g sen ( 4 s- 4 n) ( I - cos^ 0sO ^ 

In a similar fashion, the local incident angles can also be expressed in terms of the 
tilting angles 6n,4n.and the incident angle, 9 . 




e' = cos 6-n COS 0 + sLn 0tv Strv 0 cos 


sin <p' = - SLTt 0 sin ( i- cos^Q 
With the connecting relations between the local angles and the surface 
slopes known (24) can be evaluated by assuming a form for P(Zx,Zy), For o Gaussian 
surface slope distribution it may be represented as 

_ Zg (32) 


or equivalently 




2 TT 


tau Qn 

0 2tn.* 


where m is the rms slope and is assigned according to Cox and Munk's slick sea data. 

Since m is usually sufficiently small for the sea and since is 

insensitive to 0 t> for a small value of 0-n , the integration with respect to 6^ given 
in (24) con be evaluated by the method of steepest decent; if so, the result may be 
expressed as 

^ 2 >:r L *0-,= fcxn (w) 6 * -Ian (m\ 

. ~ « -n 


11.4 Modified Fresnel Reflection Coefficients 


As mentioned in section I 1.2, the Fresnel reflection coefficients should be 
modified to account for the presence of the small irregularities [12] [8j. The method 
for computing these coefficients has been discussed by Rice for horizontally polarized 
waves. Following Rice's approach also for the vertically polarized wave, we obtain 
the modified Fresnel reflection coefficients os follows 


< Rh> =Rh(0) I - 2 k COS0J I 

-Loo OQ 

clu dvr J 


-ksvnO.v) 


OO OO 

I w(U- ksunG.vj 


i - (C-bUu.+bc>- 

•Sc-n © , ■ 5— — r 

kCu.* + v*+ be) . 




^oo ■'-OO 


Cc-b) [v % be - ^ - (a + v^) sin^o) + 2 k^/7r- sLn *0 uj j d u d v 
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(35b) 


where 


Rk(0) = 

RvCG) = 


cos 0 - -si-n^Q 

COS0-+ -scTi^o 

CrCOSQ - s/gr -Si-tt^Q) 

CrCOSO +y €r - sLrv.^0 


I 


cos 0 = 7 ^ 0 + C03 0 cos 6s 



( 36 ) 


^r~ complex relative dielectric constant 


For the special case where 

Alfj* expf" — (w-*+i'’*“2kusc'n0 + k 
u- k sLn.®, v) = -4 


W( 


implifies to 1 

I<R^>|"= |Rjl'[l- 4 l<' 5 ;'c 05 e exp(-TSi>^*®)j 


|<R >1 simplifi es to 

* l<R^>l'=li 

11.5 Backscattering cross-sections 


or ^r (38) 


Substituting 0^= 6 and <#> 5 =^'ir into (36) we obtain 

RkCo) = - R^Co)= 


and 


• + /e7 

2 _ I i-yg; 


I(RkCo)>i"=I<r^co)>| =|-;^| (1-4.0 kvy) 


From (17), (19) and (3a), it follows that 


2 tan 6 
COs“-0 ^ 


(39) 


( 40 ) 


Equation (40) shows that the backscattering cross-section of the large 
undulations Is polarization Independent. 
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To find note the following two points: 

(I) Substituting 03=9 and into (30) we obtain 

COS 6s = cos 0 cos cos 6 ■+ si.n 8n cos sin 6 
Sin = - sin<f} = sinC<i>'+'^) 


(41) 


or 


as = 0' 

4*s = </> ' + 'Tr 


(2) Since 0j=0 and 4 > 5 =- 'ft' imply 63=6 and + 'ftr , from (34) 
it follows that 

'TT 

mem) ^ 0^=- birr (Tn) J 

6 s= 8 ; (42) 

<^j= 4> V-rr 

Thus from (3b) and (42), 




4*S^'!T 


Jt 


■i 


8„- tarrctn) 

eir e' 
4>+'7r 




+ o- 6 j,(e) 


0„:tarr O") 


0.= - tcm CTn) 




(43) 


H = K or V 


where 


= COS0 • 6s 


4*' + 'Tr 
65 = e' 

. 4 


(44) 


^&-lrl(6 J= COS 0 ' • ^ 6 8$> 'f^s ) 


4> j = <t> + tt 

0's = 0' 


= 4k"< f cos"0 


(€r-l) (c^r'Os'^'n^0'-+ Cr] 


(crCOse'+ye^'sLn^e' ) 




(^5) 

The complete backscattcring cross-section is, of course, given by (3) that is the sum 
of (40) and (43), 
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INCIDENT ANGLE /9 (DEGREES) 


Figure 4, Comparison of fhe angular variation of sin 0 and 35.3 exf> 9) 





HI. SELECTION OF PARAMETERS 


The surface parameters appear in the above theory are the rms slope of the 
large structures m, the standard deviation of the small structures and the 

correlation length of the small structures This scattering model can be adapted 

to predict sea returns by noting that the rms slope can be based on measurements by 
Cox and Munk [7] and that the assumed Gaussian spectrum for the small irregularities 
can approximate the high frequency po>‘tion of the sea spectrum BK ^ where the Bragg 
scatter condition applies, i.e., K = 2k sin ©. In view of the requirements of the 
composite surface theory it is reasonable to choose the oil slick sea measurements by 
Cox and Munk, since the small irregularities have been suppressed. The value of 
m is thus assigned. The value of is assigned so that the correct angular 
behavior of the Gaussian spectrum approximated BK ^ well over the angular range, 
30°< © 5 70°, i.e., BK ^ is approximated by ^ ) 

where K = 2k sin © , the Bragg scatter coridition. This is achieved by noting that 
when k/ = 2 similar behavior is realized (see Figure 4). The factor 35.3 appears 
in the Gaussian approximation to bring the levels Into agreement at ©= 60 degrees. 

The value of B must yet be assigned. 

-3 

Oceanographic Investigations indicate values of B in the range from 4.6 x 10 
-2 

to 3.26 X 10 [7, 18,19], This implies that kcr^ should lie in the range from 0.084 

to 0.24 when BK ^ is equated to ^ /n ) at ©= 60 degrees. These 

values of are consistent with the assumptions of the small perturbation theory, 
an encouraging result. The recent reports by Sutherland [11] and Pierson [20] Indicate 
that B Is a function of the wind speed. Thus, the surface parameter m.ust also be 
a function of the wind speed. It is noted that the horizontally polarized emission 
characteristic for nadir angles from zero to thirty degrees Is very sensitive to 
and hence the parameter kcr^ can be estimated by fitting the predicted emission 
characteristics to the measured data.* 

With the surface parameters established in the manner described above, both 
the emission and the backscatter characteristics may be computed and compared with 
reported measurements. 

*it appears that the wind sensitivity of B may be assigned by this technique. 
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IV. COMPARISON WITH EXPERIMENTS 


The parameter k is estimated from horizontally polarized emission 
characteristics at 8.36 GHz associated with two distinct wind speeds. The 
emission characteristics are based on an average of several of Hollinger's 
experiment runs* under similar wind stress conditions.** The vertically polarized 
emission characteristic is then computed from the estimated k <r^ . These results 
are shown in the graphs of Figures 5 and 6. The dielectric constant is based on 
data reported by Saxton and Lane [10], 

Comparison of the predictions of this emission model Indicates a 
significantly improved agreement over that predicted by a single surface model. 

Better level and trend agreement is evident for both horizontally and vertically 
polarized emissions. Sensitivity to wind speed is evident at nadir which is not 
noted with the single surface model. The sensitivity at nadir is carried by the 
modified Fresnel coefficient (see Equation (38)). 

The adequacy of the composite surface theory is further demonstrated when 
the predicted backscatter characteristics are compared with measured characteristics . 
Data at 8.91 GHz reported by Daley, et al . , *** [22] at similar wind stress 
conditions were chosen as a basis for comparison. The dielectric constant is changed 
to reflect the influence of the slightly different frequency. The comparison of 
predicted and measured' characteristics are shown in Figure 7 through 10. These 
results indicate reasonabl e agreement with measurements and improved agreement over 
the predictions of the simple geometric optics approach (single surface model). It 
is noted that the best agreement with measurements occurs primarily at larger angles 
and for vertical polarization. There is some uncertainty in the accuracy of the 
measurements near nadir [22] so lack of agreement may be anticipated there. The 
discrepancy at large angles for horizontally polarized cross sections may be attribu- 
table to receiver noise at these small cross sections. This statement is, however, 
speculative . 

*The authors are indebted to Dr. J. P. Hoi linger of NRL for making his radiometric 
measurements available to us. His data with foam and atmospheric effects removed 
are appropriate to compare with the authors calculations. 

**The authors are also grateful to Dr. V. J. Cardone of NYU for interpreting 
Hoi I inger's wind speed measurements under comparable wind stress conditions. 

***The authors are indebted to Mr. N. W. Guinard of NRL for making these backscatter 
data available to us. 
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The level of NRL data which are based on the statistical median had to be 
raised by 6 dB to realize the agreement. Valenzuela [23] indicated that the 
average cross section was about 4.6 dB above the median based on exponential 
statistics assumed for the returns. As a consequence, 1 .4 dB remains unaccounted 
for. Perhaps the remaining 1 .4 dB may be partially associated with the biases 
disclosed by Claassen and Fung. [24] 


21 


V. CONCLUSIONS 


A bIstatic two scale non-coherent scattering theory extended from 
Symyonov's paper [8] has been developed to yield the expressions for the differential 
scattering coefficients. The emission and the backscattering characteristics are then 
derived from the differential scattering coefficients in the standard way. The theory 
assumed Gaussian surface height distributions and Gaussian correlation functions for 
both scales of roughness. 

The emission and the scattering characteristics are shown to be dependent on 
the rms slope of the large undulations m, the standard deviation of the small irregularities 
, and the correlation length of the small irregularities ^ . The wind dependence 
of the first two parameters is associated with m through slick sea measurements by Cox 
and Munk, and the through the high frequency sea spectrum. The parameter^ 
is associated with the shape of the high frequency sea spectrum and can be reasonably 
chosen by fitting the sea spectrum BK ^ to the assumed Gaussian spectrum. 

It is noted that the emission characteristic for horizontal polarization is a 
sensitive measure of cr^ . Thus, cr^ is established by fitting the emission characteristic 
to measured data for different wind speeds. The parameters chosen in this way are 
then used to compute the vertically polarized emission characteristic. Good 
agreement with measured data and better agreement than a simple swface model are 
demonstrated . 

The same set of surface parameters at each wind speed is then used to compute 
the backscatter characteristics. The results except for level are shown to agree 
reasonably over all angles with NRL backscatter data under similar wind conditions. 
Comparison of these characteristics with a single parameter surface model demonstrated 
better results. 

These findings have proven that the validity of scattering theories is better 
demonstrated when both the predicted backscatter and the emission characteristics 
are compared with measurements. They have further shown that the measured emission 
and scattering characteristics with the aid of a reasonable composite surface theory 
may aid the oceanographer In Identifying the wind dependence of the sea spectrum. 
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Figure 7. Comparison of computed and measured backscattering cross-section. 
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BACKSCAHERING CROSS SECTION 



Figure 8. Comparison of computed and measured backscattering cross-section 
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Figure 10, Comparison of computed and measured bockscattering cross-sect 



REFERENCES 


[1] Nordberg, W., J. Conaway and P. Thaddeus, "Microwave Observations of Sea 

State From Aircraft," Qu arterl^y Journal of the Royal Meteorological 
Society , vol . 95, pp. ”408-413, 1969. 

[2] Nordberg, W., J. Conaway, D. B. Ross, T. Wilbert, "Measurement of Micro- 

wave Emission from a Foam Covered Wind Driven Sea," Goddard Space 
Flight Center, X-650-70-384, October 1970 (Submitted to Journal of 
Atmospheric Sciences). 

[3] Ross, D. B., V. J. Cardone, J. W. Conaway, Jr., "Laser and Microwave 

Observation of Sea-Surface Condition for Fetch Limited 17-25 m/s 
Winds," IEEE Transactions, vol. GE-8, no. 4, pp. 326-336, October 
1970. 

[4] Hollinger, J. P., "Passive Microwave Measurements of the Sea Surface, " Journal 

of Geophysical Research, vol. 75, no. 27, pp. 5209-5213, September 
1970. 

[5] Hollinger, J. P., "Passive Microwave Measurements of Sea Surface Roughness," 

IEEE Transactions, vol. GE 9, no. 3, pp. 165-169, July 1971. 

[6] Stogryn, A., "The Apparent Temperature of the Sea at Microwave Frequencies," 

IEEE Transactions on An tennas and Propagation, vol. AP-15, no. 2, 
pp. 278-286, March 1967. 

[7] Cox, C. and W. Munk, "Statistics of the Sea Surface Derived from Sun Glitter," 

Journal of Marine Research, vol. 13, no. 2, pp. 198-227, February 1954. 

[8] Semyonov, B., "Approximate Computation of Scattering of Electromagnetic 

Waves by Rough Surface Contours," Radio Engineering and Electronic 
Phys., vol. 11, pp. 1179-1187, 1966. 

[9] Hollinger, J. P., Private Communication. 

[10] Saxton, J. A. and J. A. Lane, "Electrical Properties of Sea Water," Wireless 

Engr . , vol. 29, p. 269, October 1952. 

[11] Sutherland, A. J., "Spectral Measurements and Growth Rates of Wind-Generated . 

Water Waves," Stanford University, Dept, of Civil Engineering, Technical 
Report 84, August 1967. 

[12] Peake, W. H., "Interaction of Electromagnetic Waves with Some Natural 

Surfaces," IRE Transaction, vol. AP-7, Special Supplement, pp. S324- 
S329, December 1959. 


29 


[13] Rice, D. O., "Reflection of Electromagnetic Waves from Slightly Rough 

Surfaces," Communications in Pure and Applied Mathematics," vol. 4, 
pp. 361-378, February 1951. 

[14] Fung, A. K., "Theory of Cross-Polarized Power Returned From a Random 

Surface," Appl . Sci. Res. 18, pp. 50-60, August 1967. 

[15] Stogryn, A., "Electromagnetic Scattering From Rough, Finitely Conducting 

Surfaces," Radio Science , vol. 2, (New Series), no. 4, pp. 415-428. 

[16] Valenzuela, G. R., "Depolarization of EM Waves by Slightly Rough Surfaces," 

IEEE Transaction s on Antennas an d Propagati on, vol. AP-15, no. 4, 
pp. 552-557, July 1967. 

[17] Barrick, D. E., and W. FI. Peake, "Scattering from Surfaces with Different 

Roughness Scales: Analysis and Interpretation," Res. Rep. BAT-197A- 
10-3, Battel le Memorial Institute, Columbus Laboratories, Nov. 1, 1967. 

[18] Pierson, W. J., Jr., "A Proposed Vector Wave Number Spectrum for a Study 

of Radar Sea Return, Microwave Observations of the Sea," NASA/NAVY 
Review, SP-152, pp. 251-282. 

[19] Phillips, O. M., "The Dynamics of the Upper Ocean," Cambridge University 

Press, London, p. 120, 1966. 

[20] Pierson, W. J., Jr., Private Communication. 

[21] Cardone, V. J., "Specification of the Wind Distribution in the Marine Boundary 

Layer for Wave Forecasting," Ph. D. Thesis, New York University, 
Geophysical Science Laboratory, Technical Report 69-1, December 1969. 

[22] Daley, J. C., J. T. Ransome, Jr., and J. A. Burkett, "Radar Sea Return — 

JOSS I," Naval Research Laboratory, Report No. 7268, May 1971. 

[23] Valenzuela, G. R., M. B. Laina and J. C. Daley, "Ocean Spectra for Fligh 

Frequency Waves from Airborne Radar Measurements," Journal of Marine 
Research, vol. 29, no. 2, May 1971. 

[24] Claassen, J. P. and FI. S. Fung, "The Wind Response of Radar Sea Returns and 

Its Implication on Wave Spectral Growth," University of Kansas Center 
for Research, Inc., Technical Report 186-5, September 1971. (Available 
upon Request) 


30 


NASA-Langley, 1973 20 CR-2326 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE $300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
45! 



POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


'*The aeronautical and space activities oj the United States shall be 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, a lasting contribution to existing 
knowledge'. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Informa"tion on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from; 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

ADMINISTRATION 


NATIONAL AERONAUTICS AND SPACE 

Washington, D.C. 20546 





